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Abstract 

We study models of continuous time, symmetric, Z rf -valued random walks in random 
environments. One of our aims is to derive estimates on the decay of transition probabil- 
ities in a case where a uniform ellipticity assumption is absent. We consider the case of 
independent conductances with a polynomial tail near and obtain precise asymptotics 
for the annealed return probability and convergence times for the random walk confined 
to a finite box. 

1 Introduction 

We study continuous time, irreducible, symmetric, nearest neighbor random walks in random 
environments on Z d . Our aim is to derive estimates on the decay of transition probabilities in 
the absence of a uniform ellipticity assumption. 

The paper has four sections (other than this introduction). Sections 2 and 4 deal with 
the decay of the mean or annealed return probability. In Section 2, we consider quite general 
reversible random walks in a random environment and we establish a comparison lemma for 
the annealed return probability. The proof is based on a trace formula (in fact an extension 
of the trace formula for central probability for random walks on amenable groups, see [9]). In 
Section 4, we derive sharp bounds on the decay of the annealed return probability from direct 
investigation of traces and eigenvalues when the rates are i.i.d. random variables chosen from a 
law with polynomial tail near 0. We then prove that one might get the classical t~ d l 2 decay or 
a slower decay of the form t -7 , where 7 < d/2 is related to the tail of the law of the rates near 
0. In Section 5 we deal with the quenched decay and obtain a partial result (Theorem 5.1) that 
nonetheless establishes a difference with respect to the annealed decay for small values of 7. 

In Section 3, we discuss finite volume random walks taking their values in a torus. We 
obtain some quenched estimates on convergence times when the random rates are i.i.d., chosen 
from a law with polynomial tail near 0. These follow from sharp bounds on the spectral gap. 
In particular we prove a universal lower bound for the spectral gap of a symmetric random 
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walk on a torus of side length N (Proposition 3.20 below) which allows to separate the effects 
of the usual diffusive N~ 2 factor and the contribution of small values of the rates. 

The paper is written in such a way to ease independent reading of the different parts at 
the cost of some repetition. Sections 2 and 3 are self-contained; only the spectral gap lower 
bound (3.35) from Section 3 is needed to proceed through Section 4. 

2 A comparison lemma for the annealed return proba- 
bility 

We study a family of symmetric, irreducible, nearest neighbors Markov chains taking their 
values in Z d and constructed in the following way. Let Q be the set of functions w : Z^xZ^ -> M + 
such that u(x,y) > iff x ~ y, and u(x,y) = cu(y,x). (y ~ x means that x and y are nearest 
neighbors.) We call elements of Q environments. 
Define the Markov generator 

^/(x)=X)w(x,i/)[/(y)-/(x)]. (2.1) 

y~x 

As usual, {X t , t G P+} will be the coordinate process on path space (Z d ) R+ and we use the 
notation P£ to denote the unique probability measure on path space under which {X t , t G M + } 
is the Markov process generated by (2.1) and satisfying X = x. Under P£, X = x; then 
the process waits for an exponentially distributed random time of parameter Yly^x u(x, y) 
and jumps to point x\ with probability <jj{x,x\) / ^ x uj{x,y)\ this procedure is then iterated 
choosing independent hoping times. Equivalently, one can define P£ using the theory of sym- 
metric Dirichlet forms, see [4]. The reference space is then L 2 (Z d ), equipped with the counting 
measure. For functions / and g with finite support, let 

V*(f,9) = \ E <*frv)\f(*)-f(v)M*)-9(v)]- 

The bilinear form T> u is closable and its closure is a regular, symmetric Dirichlet form. Thus, 
there exists a Hunt process associated to V w . Note that points have non zero capacity. There- 
fore, the measure P£ is uniquely determined by V w . It is easy to prove that both constructions 
yield the same law P£. 

Since w(x,y) > for all neighboring pairs (x,y), X t is irreducible under P£ for all x. The 
counting measure on Z d is reversible because we have assumed that u(x,y) = cu(y,x). 

We now choose the rates oo at random, according to a translation invariant law Q on VL. 

In the sequel Q.IP£ will be used as a short hand notation for the annealed law defined by 
Q.P£[.] = J P£[.]dQ(u). Note that X t is Markov under P£ for any cu, but is not Markov 
anymore under Q.P£ for nontrivial Q. Let P w = P£ and Q.P W = Q.P£. 

We are interested in estimating the decay of the annealed return probability Q.P^ [X t = 0] , 
as t tends to +oo. 

As a subset of (R + ) zdxZd , Q is a partially ordered set. By duality, one can define a partial 
order on the set of probabilities on f2 in the following way. Given two probabilities, Q and Q', 
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we say that Q' > Q if, for any measurable, bounded, increasing function / : Q — > R, we have 
Q'(/) — Q(/)- (/ is increasing if, whenever uj, uj' e Q satisfy uj'(x, y) > uj{x, y) for all x, y, then 

iV) > /HO 

Remark 2.1 T/ie function uj — » P^L^Q = 0] noi monotonous in uj. It is clearly not increas- 
ing. It is also not very difficult to find subgraphs ofL d for which the removal of an edge decreases 
the value o/P^L^Q = 0] (left as an exercice), which implies that the function uj — > P w [X t = 0] 
is not decreasing. 

Lemma 2.2 Let Q and Q' be two probabilities on Q such that Q' > Q. Assume that for 
Q' + Q-almost all uj, the Markov chain X t is conservative under P w . Then, for all time t, we 
have 

Q'.P^Xt = 0] < Q.^[X t = 0]. 

Proof We prove that Q.P w [X t = 0] can be written as a supremum of the Q-expectation of 
decreasing in uj functions. More precisely, let B N = [-N, N] d be the box centered at the origin 
and of radius N. Let Q W,N be the restriction of the operator Q w to Bn with Dirichlet boundary 
conditions outside (that is, Q^' N is the generator of the process which coincides with the 
one given by Q w until the latter process leaves B^ for the first time, and then it is killed). 
Then — Q W,N is a positive symmetric operator. Let {/if(B N ),i e [1,#B N ]} be the set of its 
eigenvalues labeled in increasing order. We shall prove that 



.W»[X t = 0] = sup- 



n #B N 
Let 



E 



(2.2) 



^' N U\g) = l E <*(*>v)\f(x)-f(y)Mx)-9(y)]+ E faM*) E "( x >v) 

x,yeB N " yt B N 

be the Dirichlet form of — Q U ' N . From the min-max caracterization of ^(Bn), we have 

W t T> \ ^ (/' /) 

Hi (B N ) = max mm — — — , 

where the 'max' is computed on choices of % — 1 functions defined on B^ and the 'min' is 
computed on functions / such that, for all j G {1, . . . ,i — 1}, Yl x eB N f( x )fj(y) = 0- For any 
function /, £ w,N (f,f) is clearly increasing in uj, therefore for given N, and i, ^(Bn) is an 
increasing function of oj and ^.e^^'* is decreasing in uj. Thus (2.2) implies the lemma. □ 

Proof of (2.2) Let r N be the exit time of X t outside B N . Note that £\ e - ""^)* is just the 
trace of the semi-group of the process X t killed when leaving the box B N , i.e., 

Y^ e -MB N )t = J- FZ[X t = x;t<r N ]. 
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We compute Q.Pq [X t = 0] using that, from the translation invariance of the probability Q, 
we know that Q.F^[X t = x] does not depend on x. Therefore 

q.p£[x, = o] = ^i- Q- W = A > E ^- p * ^ = x ' t< ^ 



xeB N 



N 



-H?(B N )t 



proves the lower bound. 

As far as the upper bound is now concerned, note that 

Q.vS[x t = o] = -L_ Q-^I** = *] 



= 4J~R~ E Q-^t** = ^ * < TN ^ + ^777- E ^ = ^ * ^ ^ 



xeBM 



x&B N+k 



xeB N 



We have 



E' 



-^f(BN+k)t 



< #B N+k sup 



1 



M #-B 



a; 



(B M )t 



Let ra t be the number of jumps the process X t performs by time t. For x G -Bat, under P£, 
* > tjv+a: implies that n t > k. Therefore 

E Q-^[t>T N+ k] < E Q-^x[nt>k} 

x£Bn x6-Bjv 

= #S JV Q.P w [n t >fc], 

using the translation invariance in the last equality. 
So far, we have obtained the bound 



E< 



-H?(B M )t 



#B N m #B M 
First let N tend to +oo, then let k tend to +oo to deduce that 

1 



+ Q.P a >t > k]. 



.W%[X t = 0] <sup 



M #B M 



E 



DJ*" = +oo 



Now the conservativeness assumption and the fact that there are no instantaneous points 
of X t in Z d imply that P w [n t = +oo] = Q-a.s. □ 



4 



3 Times of convergence to equilibrium of random walks 
on the torus 

Let S N be the discrete, <i-dimensional torus of side length N. When convenient, we consider 
Sn as a subset of Z d . We construct a family of Markov chains taking their values in Sn- Let 
uj : Sn — » R+, and define the Markov generator 

C"> N f(x) = 5>(s) A u(y)] [f(y) - f(x)}, (3.1) 

where the sum is over sites y which are nearest neighbors to x (relation that is denoted y ~ x). 
Let {X t , t G R + } be the process with distribution P^'^ generated by (3.1) and the condition 
Xq = x. Since uo(x) > for all x, X t is ergodic under P^-^ for all x. The unique invariant 
probability measure is the uniform law, denoted by i]n- Furthermore, 77 at is reversible. 
We choose the family {u(x), x G Z d } i.i.d. according to a law Q on (M* + ) zd such that 

uj(x) < 1 for all x; (3.2) 
Q(w(0) < a) ~ a 7 as a | 0, (3.3) 

where 7 > is a parameter. 

Remark 3.1 VKe note t/ioi this generator has the same form as Q w in (2.1) by making u(x, y) = 
uj(x) A uj(y), but for a process in finite volume. We could have defined uj on edges, instead of 
points, as in the previous section, with i.i.d. values for different edges, and the same technique 
would apply, with similar results, and heavier computation. 

Remark 3.2 lfu(0) were a Bernoulli random variable, then we would have a random walk on 
a (independent, site) percolation cluster (provided we started in an occupied site). See [6]. 

Our main results refer to the following convergence times. For e G (0, 1), let 

T?' N = ini{t > : sup sup K' N [f(X t )} - rj N (f)\ < e} (3.4) 

xes N |/|<i 

T^ N = mi{t > : sup sup \V%f\f(X )g(X t )] - Mf)Mg)\ < e}, (3.5) 

l/l<i|sl<i 

where E%' N is the expectation with respect to P^-^ and E^(-) = jE^ N (-) dr] N (x). 

Remark 3.3 The first convergence time is a worst-case one, that is, it is the longest conver- 
gence time among all initial conditions. The second one is an average convergence time among 
all initial conditions (under uniform weighting). 

Remark 3.4 Clearly, T% ,N < T± ,N for all uj. 

Theorem 3.5 For all 7 > and e G (0, 1), we have Q-a.s. 

limsup , \ T < 2V-, 3.6 
^oo P logiV " 7 

lim . nf fogrr^ 2y d_ 

N^oo fog N 7 
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Theorem 3.6 For all 7 > and e G (0, 1/4), we have Q-a.s. 

limsup ° g 2 < 2, (3.8) 
iv-»oo \ogN ~ y ' 

liminf °° 2 > 2. (3.9) 

N^oo \ogN 

In fact, for all e G (0, 1/4), there exists a constant c > such that for all uj 

liminf N~ 2 T^ N > c. (3.10) 

Remark 3.7 // uj{0) were bounded away from zero, that is, if uj(0) > c\ Q-a.s. for some 
constant C\ > 0, then limsup^^^ N~ 2 T± ,N < c 2 Q-a.s. for some constant c 2 > 0. 

Remark 3.8 Theorems 3.5 and 3.6 establish that Q-a.s. 

lim i^S = 2 V * and lim ^ = 2. (3.11) 

A^oo logiV 7 iV^oc logiV 

We £/ius nave distinct asymptotic behaviors of T^' N and T% ,N when d/j > 2. ^4 heuristic 
argument to justify that follows. When d/7 > 2, T^ ,iV ; as a worsi case convergence time, is 
greater than or equal to the convergence time starting at a site with minimal uj, whose order 
is clearly smaller than or equal to N~ d ^ . On the other hand, choosing a site uniformly at 
random as a starting point will miss the low uj sites and, starting at high uj, the walk will get to 
equilibrium faster than it will get to any low uj site. It will be as if there were no low uj sites, 
and that means T2 is of order N 2 (see Remark 3.7). 

From now on, we shall drop the W in some of our notation. For example, we use the short 
hand notation S = Sjy. 



3.1 Proof of (3.10) 

Let A = {x = (xi, . . . ,x d ) G S : x l G [0,N/2]}, T A = M{t > : X t G A} and, for A > 0, 
h%(\) = Ef N (e- XTA ). Choosing / = 1 A and g = l A c, we have 

sup \E»f[f(X )g(X t )\ - VN (f) VN (g)\ > Vn (A) Vn (A c ) - P£f (X £ A,X t G A) 
l/l,lffl<i 

> rj N (A)rj N (A c ) - P^(X i A,T A < t) > rj N (A)rj N (A c ) - inf m (l A ch" (A))e A *. (3.12) 

We now estimate r)N(l-A c h w (\)). We will compare with the case uj = 1, which corresponds to 
the usual random walk on Sjy. The Dirichlet form of X t is given by 

£"> N (f, /) = ^E M 1 ) a <*(vWW - f(y)) 2 - ( 3 - 13 ) 

x~yes 

It is clear that £ u '' N (f, f) is nondecreasing in (the natural partial ordering of) uj. We have also 
that, for A > 0, 

\r] N (h"(\)) = inf £"> N (f, f) + \ VN (f 2 ). (3.14) 
/U=i 
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Since £ U ' N (/, /) < S 1 (f, f), where 1 is the identically 1 vector indexed by S, we have that 

VNih^iX)) < r] N (h 1 (X)). (3.15) 

Since T A is a hitting time for an ordinary rate 1 random walk on Z under P 1 , the invariance 
principle yields that for all A > 

^(/ i 1 (iV- 2 A))^i + 0(A) (3.16) 

as iV — > oo, where 0(A) — > as A — > oo. We also have that rj N (h^{X)) = r] N (A) + ^(l^c/i^A)) 
and ^(A) -> 1/2 when N ^ oo. Thus, from (3.14), (3.15) and (3.16), 

limsup7^(/^(iV- 2 A)) = \ + limsup?7 A r(lAc/i a; (^ _2 A)) < J + 0(A) (3.17) 

AT— >oo 2 AT— >oo 2 

and it follows that 

W (l Ac ^(iV- 2 A)) < 0(A). (3.18) 

We conclude that 

liminf sup \E^ N N [f(X )g(X cN2 )]-r ]N (f)r ]N (g)\>]-e ( l>(l/c). (3.19) 
l/l,lsl<i 4 

Since 0(l/c) — > as c — > 0, we get that for all e < 1/4, lim inf at^qo N~ 2 T%' N > c*, where c* is 
any positive constant satisfying 0(l/c*) < (1/4 — e)/e. □ 

3.2 Proof of (3.8) 

We make use of generalized Poincare inequalities [7], which we recall now. Let B denote the 
set of nearest neighbor bonds of S, i.e., B = {(x, y) : x,y G S,x ~ y}. For x,y & S, define 
r"(b) = N- d (u(x) A u(y)), if b G B, and r M (b) = 0, otherwise. The Dirichlet form of £"> N on 
L2(S,t]n) can be written as 



2 

fees 



where (4/ = /(x) — /(?/) and the sum ranges over b — (x,y), x,y G 5. 
For p G (0, 2), let g be such that 1 + 1/g = 2/p and 

r^ N {p) = sup ^ (/2)2/P 2/ . (3.20) 

m, mU )=o £^ N (f,f)\\f\& q 

We then have 

T^ N < qe-Vi^ip) (3.21) 

for all p G (0,2). 

Remark 3.9 In the notation of [7], T^^ip), as defined in (3.20), equals l/IC UJ (p). 
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For all x,y G S, let n Xjy be a nearest neighbor path from x to y and let t = sup xy \n x ,y \ be 
the length of the longest path. 

Consider now a partitioning of S = B U G and let 

= T. b ={ x ,y)^G( d bf?VN{x)T] N {y) 

Lemma 3.10 



t"> n ( P ) < 2 2/q 3 2/p r] N (B) 2/p t sup — - + 2 2 '\ 

b=(x,y):x~y r \U) 



' G 



In the next lemmas, for given < £ < 1, we choose G as the largest connected component 
of the set {x : u(x) > £} (following a deterministic order in case of ties). 

Lemma 3.11 For £ > small enough, there exists a positive number c\ that depends only on 
d such that Q-a.s. 

liminf N 2 /t%' n > fci. 

JV— »oo 

Lemma 3.12 For £ > 0, there exists a number c 2 (£) sncn inat c 2 (£) — > as £ — > and Q-a.s. 

lim sup ?7jv(-B) < c 2 (0, 

wnere B = S\G. 

Lemma 3.13 There exists a finite number c 3 depending only on d such that Q-a.s., for all N 
large enough 

in£u(x) > N~ C3 . 

In the proof below, we will see that c 3 can be taken as ^ + e for arbitrary e > 0. 
We postpone the proofs of the above lemmas until after the proof of (3.8). 

Proof of (3.8). With uj and £ > fixed, we choose N big enough so that the conclusions of 
Lemmas 3.10, 3.11, 3.12 and 3.13 hold. Then, using also (3.21), 

T?' N < qe~ 1/q T"' N (p) < ^{[12e- 1 c 2 (0] 2/p iV 1+d+Cs + (ciO -1 ^- 1 ) 2 ^ 2 ]}. (3.22) 

Assuming that £ is small enough, let p satisfy 

2 = (rf-l + c 3 )logiV + log(ciO 
P logc 2 (0 _1 - log 3 

With this choice, the two summands in the expression within braces in (3.22) are equal and 
thus (3.22) equals 

gffc^Wexp lQg(Cl ° cxp { U + (^- 1 + ^)log(4e- 1 ) l 1 

2 (Cie) 6XP logMO)- - log3 6XP \ L 2 + log^))" 1 - log3 J l0gN S ■ (3 - 23) 



Combining (3.22-3.23), we get 

lim sup 



logT?'" n , (rf-l + c 3 )log(4e- 1 ; 



< 2 + 



logiV log(c 2 (0)- 1 -log3 • 

Since this holds for all £ > sufficiently small and c 2 (£) — > as £ — > 0, the result follows. □ 

Proof of Lemma 3.10. This is very similar to the results of part III in [8]. We estimate the 
three terms in the decomposition 



Mf) 



( \ E + E + \ E ) (/(*) " f(v)) 2 m(x)m(y) =-.i + n + m (3.24) 

V x,ydG xeG,yeB x,yeB/ 



in turn. 



7 < 



Holder 



x,y&G 



p/2 



< 2^||/||^ £ (/(*) - /(y))^(x)^(y) 



\x,yeG 



< 2^11/11^ (2^(/,/)/r^) 



p/2 



(3.25) 



// < (211/IU) 2 ^ £ (/(*) - f(y)) p VN{x) VN {y) 

x&G,y&B 



Holder 



< (2||/|kr"| £ (/(*) - /(y) W^dl) (^(G)^(S)) 1 " 

\x£G,y&B J 

< (w\u 2 - p ( E IE 
/ V 



< (211/iu 



2-p 



E E M x )Vn{v) 



b x<EG,y(£B: 



MB) 



i-p 



Holder , 

< (2||/||oo) 2 - p (£ V(6) 



p/2 



/ 




/ 


( tv /2 




Mz)Mv)\ MB) 1 '* 




x£G,yeB: J J 



if \ p/2 
< (2||/|| 00 ) 2 ^(2^(/ ! /)) p/2 ^sup(r-(6))- 1 J m {B) 



p/2 



X 



V 



/ \ P/2 

< ^- p/2 \\f\\^ P [£ u '' N (f,f)sMr^b))- 1 j (tvN(B)f 2 VN (BY-^ 

( \ p/2 
= 2 1 -^||/||L-V(5)^' JV (/,/)sup(r-(6))- 1 rj , (3.26) 

where the last inequality follows from 



Vn{x)vn{v) = Y \ n *,y\VN{x)VN{y) < trj N (B). 

Similarly, 



b xea, y eB: x£G,yeB 



III < 2 1 - p / 2 \\f\\ 2 ^r ]N (B) ^(/,/)sup(r-(6))- 1 rj . (3.27) 
We conclude from (3.24), (3.25), (3.26) and (3.27) that 

Thus, 

r^ N (p) < U2- p / 2 r ]N (B)\sup(r"(b))- 1 tj + 2^ (t£ n J \ 

< 2I- 1 |(3^(5)) 2 / p sup(r w (6))- 1 r +7# Jv j 

= 2l +1 |(3r ?iV ( J B)) 2 ^sup(r-(6))- 1 r + r^|. □ 

Proof of Lemma 3.11. Since w(-) > ^ on G, we have 
where 

1 _ Y,b=(x,y)eGxG( d l>f) 2 (# G )~~ 

T ° : " 5S £ 6eGxG (4/) 2 (#G9-i 

is the inverse of the spectral gap for the ordinary rate 1 random walk on G. From Cheeger's 
inequality, we get that 

r G < 8E 2 G , 

and therefore 

r% N < Sr 1 ^, (3.28) 
where the isoperimetric constant E G is defined by: 

#A#G\A 
:= sup 7777-77^ . , 
acg #G#d G A 
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where d G A = {(x, y) : x ~ y, x G A, y G G \ A} is the bond boundary of A with respect to 
G. The statement of the Lemma will thus follow if we can prove that NEq 1 is bounded from 
below for large N by some constant that only depends on the dimension. We shall rather show 
that 

^Q(~ G >aiV)<oo, (3.29) 

N 

for some a. One then uses the Borel-Cantelli Lemma to deduce from (3.29) that, Q.a.s., for 
large N, we have E G < aN and therefore, as follows from (3.28), r G ' N < 8^a 2 N 2 . 

Following [6], Subsection 3.1, we note that we can restrict ourselves to connected A's such 
that G \ A is connected. 

Since #d G A > 1, we have f§jf^ < %N as soon as #A < |7V or #G \A< f N. Thus we 
may also assume that #A > f N and #G \ A > f TV. 

The same argument as in [6], Subsection 3.1, based on the classical isoperimetric inequality 
on S, shows that (3.29) follows from 

VQ(sup— 77 : — —>a]<oo. (3.30) 

In (3.30), B = {(x, y) : x,y G S,x ~ y} denotes the set of nearest neighbor bonds of S. The 
sup is computed on *-connected sets F C B such that #F > a\N~, for some constant ct\ 
that depends on a and the dimension. 

Given such an F, choose a subset, say F, such that b = (x, y) ^ b' — (x', y') G F =>• x ^ x' 
and y ^ y'. Since any point has at most 2d neighbors and #F > aiN~ , we may assume that 
F > «2#F, for some positive a 2 - 

Now, for all A > 

Q(#{(x,y) G F; a;(x) > ^(y) > £} < #F/a) 

< Q(#{(z, y) e F; cu(:r) > £, wfo) > £} < #F/a) 

< e ^#^(l_ 7r 2 + e -A 7r 2 ) #F 

< e ^(l_ 7r 2 + e -A 7r 2 ) a 2 #F 



where 7r = Q(ou(x) > £). 

By the above inequality, and the fact that the number of distinct *-connected subsets F 
with #F = n is bounded above by N d e a:in for some 0:3 [10], we get 

jtE. S„l <- /V d V „[a3+Aa- 1 +a 2 log(l-7r 2 +e- A 7r 2 )]n 



I sup — r ^ F — — > a ) < N d V 

n>aiiV" 



d-i 



d-l 
n>ct\N d 
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where 0:4 := —[a 3 + \a 1 + a 2 log(l — n 2 + e x n 2 )] > 0, provided we choose A and a such that 
a 3 + \/a < \a 2 and £ < £ , for £ close enough to 0, depending on a, A, a 2 , a 3 and 7 only. □ 

Proof of Lemma 3.12. Consider the site percolation model on Z d where a site x is occupied 
if u>(x) > g Let ^0 be positive and satisfy Q(u(x) > £0) > Pc the critical density for the 
a.s. appearance of an infinite connected component C. Then, if £ < £ , C exists a.s. Let 
C N = C fl where 

Sn is iS/v viewed as a subset of Z d (that is, without the boundary 
identification), say, Sn = (—N/2,N/2] d n'Z d . Let Cn be Cjv viewed as a subset of the torus Sn 
(that is, with the boundary identification). Then, it follows by standard ergodicity arguments 
that liniTv^oo r} N (C N ) = 6(g) := Q(0 E C) Q-a.s. Since 6(g) -> 1 as £ -> (a well known 
result [5]), the result would follow if Cn were connected, which it is not necessarily. 

Consider then Cn '■= Cat-l^j- We claim that CV is connected in Sn, and thus also in 

Sn, for all large enough Q-a.s. Indeed, in the event that Cn is not connected in S N , there 
exist two sites at the boundary of S N _^^ N -^ that are connected to the boundary of Sn but are 

not connected to one another. This implies that there exists a site x at the boundary of Sn 
whose (occupied) cluster (in Sn) has a boundary (of vacant sites) of size at least [V^J- Now, 
the (bond) boundary of any finite cluster of a site in Sn can be identified with a surface of 
plaquettes around the given site, each plaquette crossing orthogonally a boundary bond. For 
each such plaquette, there corresponds thus an inner occupied site and an outer vacant one. 
For a given such surface of plaquettes of size (total number of plaquettes) n, there is at least 
n/(2d) distinct outer vacant sites (since a vacant site can not be adjacent to more than 2d 1 In 
the case of x, the surface of plaquettes will intersect the boundary of Sn in a closed curve. It 
will also have to cross the region between the boundaries of Sn and S N __ . For this reason 

it will contain at least [\/N\ plaquettes. 

From the arguments in the latter paragraph, we get the following estimate. 



/jv is not connected in Sn) 
— Z~2x£ds N Xl rarOT j^j Q( & U the sites at the outer boundary of T are vacant), (3.31) 

where the latter sum above is over surface of plaquettes T around x. The number of distinct 
such surfaces which have size n can be estimated to be exponential in n [10]. Proceeding with 
the estimation we get that the right hand side of (3.31) equals 

Ylx&ds N 2~2n>i^\ X^ rar ° u _? d:c: Q( an the sites at the outer boundary of T are vacant) 

< N d Z n >iSNi e-[Q(u;(0) < 0] n/(2d \ 

where v depends only on d. Thus, by taking < £ < £0 small enough, the probability in the 
left hand side of (3.31) can be made summable and the claim at the beginning of the previous 
paragraph follows by Borel-Cantelli. The lemma then follows. □ 

Proof of Lemma 3.13. We will prove that Q-a.s. 

log inf x u(x) _ _d 

N~*oo log iV 7 



1 Actually, n/(2d~ 1) is a better bound. 
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For that, let c < d/7. Then 

Q(inf > N- c ) = [Q(lo(x) > N~ c )f d < (1 - aN-"*)"" < e ~ ClA ^ C7 (3.32) 

for iV large enough and some constant Ci. Thus the Borel-Cantelli lemma implies the upper 
bound in (3.32). 

Now, let c > d/7. For e k < N < e k+1 , we have 



Therefore, 



inf u;(a;) > inf u;(x) A inf u(x). 
xes N xes ek x£S ek+1 \S ek 



3N G [e fc , e fe+1 ) : inf u(x) < TV" 

x£S N 



< Q inf < e' ck + Q inf cu(x) < e" 

\^s ek J \xes ek \s ek 

= (1 - (1 - Cl e-^ k f dk ) + (1 - (1 - Cie -c 7 fc)e d(fc+1) -e dfe ) < ^ e -(cr-«0* (3.33) 

and the result follows from Borel-Cantelli and the summability of the probabilities on the left 
hand sides of (3.32) and (3.33), implied by their right hand sides. □ 

3.3 Proof of (3.7) 

From T?' N > T%' N and liminf^oo N~ 2 T^' N > c Q-a.s., we deduce that liminf^oo N~ 2 T?' N > 
c Q-a.s. and, thus, lim inf jv^oo log T?' N /log N > 2 Q-a.s. 

We argue now for the inequality liminfjv^oo log T^'^/ log iV > d/7 Q-a.s. Let x G S. During 
an exponential time of parameter J2 y - y ~ x u (y) A ^{x), the process X starting at x stays still. 
Therefore, 

sup\E" x > N f(X t )-r) N (f)\ > ^ N (X t = x)-N- d 
l/l<i 



i.e., 



Therefore, 



> e -fEy~x u (.v) Aoj ( x ) _ ]V^ d > e -2 ^^)* _ N~ d , 



T?> N > ± S u P cu(xy 1 log(e + N~ d )-\ 



logT"' N logsup^rr) 1 

— 1„„ AT °V 



log N ~ logN 
Now, let < 5 < 1 be arbitrary. 



log sup u(x)' 1 < (l-5)-logV| = [Q(u(x) > N-( 1 - s ^)] Nd < [1 - N-^-^Y 4 , 

7 / 

for any 1 > 5' > 5, provided TV is large enough. Thus, the above probability is summable in N 
for any S > 0, and the result follows by Borel-Cantelli. □ 



13 



3.4 Proof of (3.6): Spectral gap estimates 

Let 




where / is any function uniformly bounded by 1. Thus 




Using a formula of Saloff-Coste (see Theorem 3.2.3 in [12]), we get 




(3.34) 



where u)(b') = cu(x') Au(y') for b' = (x 1 , y') G B, W : B — > (0, oo) is an arbitrary weight function, 
{^x,y '■ {x,y) E S x S} is an arbitrary complete set of paths (ir x ,y is a path with end points x 
and y), for an arbitrary path 7r in S, \ir\ w = ^2 bE7r 1/W(b), and Af(b, b') := y) G 5 x S : 

It remains to estimate the right hand side of (3.34). The key point here is the choices of the 
weight function and the complete set of paths. Roughly speaking, the latter will be taken in 
such a way that no path in it has interior sites with low values of u; and the former will give 
low weight to bonds with low values of uj. We are precise next. 

Definition 3.14 Given e > 0, a site x G S will be called e-good if uj{x) > A r ~ e . Otherwise, it 
will be called e-bad. A bond b = (x,y) G B will be e-good if x and y are e-good. Otherwise, it 
will be called e-bad. 

Definition 3.15 Given L > and a path n G S connecting given sites x,y, a site z in n will 
be called an L-interior site of n if \\z — x\ |oo, | \z — y\ |oo > L. 




Definition 3.16 Given L, e > and T, a set of paths of S, T will be called (L,e)-good if all 
the paths of T have all their L-interior sites, if any, e-good. 

We now construct for every iV a complete set of paths for Sn which will turn out to be 
almost surely (L, e)-good for all large enough iV and which will have other properties leading 
to the validity of (3.6). 

We start with an auxiliary set of paths. 

Definition 3.17 For x,y G S, let r\ x ^ be the path given by moving sequentially in the 1-st, 
2-nd,..., d-th coordinate direction one step at a time, along the longest segment (and according 
to an arbitrary predetermined order in case of a tie), from x to y, until the coordinates are 
successively matched. 
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For example, if d — 3, N — 100, Sn = {1, 2, . . . , 100} 3 (with the boundaries appropriately 
identified), x = (1, 1, 1) and y = (2, 20, 80), then r) XtV = 7i U72 U73 is the union of the segments 

71 = {(1,1,1) = (100, 1,1), (99, 1,1),..., (3, 1,1), (3, 1,1)} 

72 = {(2, 1,1) = (2, 100,1), (2, 99,1),..., (2, 21,1), (2, 20,1)} 

73 = {(2, 20,1), (2, 20, 2),..., (2, 20, 79), (2, 20, 80)}. 

Now for L > we define the L-sausage Sl = Sl(x,u) with base r\ x ^ and width L as follows. 
We suppose N > 3L. Let i\, i 2 , ■ ■ ■ , ik, 1 < k < d be the coordinates where x differs from y 
in increasing order, so that rj x>y is the union of the segments 71, . . . , 7fe, each of length at least 
N/2, with 7j parallel to the coordinate direction i. If k < d, then let i* = min{i : 1 < % < 
d and i ^ ik, 1 < k < d} and 

Sl = ■ ■ ■ , Zi*-\,Wi*, z,i* + i, . . . , Zd) '■ 

Zi* < Wi* < Zi* + L — 1, (Zi, . . . , Zi*-i, Zi*, Zi* + i, . . . , Zd) G 1] x ,y}- 

If k = d, then let 

s 'l = {{wi,z 2 ...,z d ):z 1 <wi<z 1 + L-l,(zi,...,z d )eU k j= f{j}, 
S 'l = {(wi, z 2 ■ ■ ■ , z d ) : zi -L+ 1 < wi < z 1 ,(z 1 ,. . . ,z d ) E U* =2 7j}. 

Now let i?i be the uniquely defined rectangle with base 71 and width L such that either Ri r\S' L 
or i?i fl S'l is a L x L square (one and only one of these possibilities occurs). In the latter case, 
S L = Ri U S'l; in the former one, S L = RiU S' L . 




Figure 1: Thick polygonal is rj Xjy ; rectangle delimited by dashed lines is R\\ strip delimited by 
thin lines is S'[ 



Remark 3.18 Notice that Sl can be seen as either a single bidimensionaP strip of length at 
least N/2 and at most dN and width L, when k < d, or the union of two such strips (one of 
which is the rectangle R\), when k — d. 

2 Even if living in fc-dimensional space. 
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Given e > and a strip S of length at least N/2 and at most dN and width L, we consider 
the site percolation model in S in which a site is open if and only if it is e-good and define 
the event A s = As(L) that there exists an open path connecting the two smaller sides of S 
(within S). Then one argues as usually that A c s is the event that there exists a *-closed path 
connecting the two larger sides of S (within S). It is clear that A${L) C As(L') if L < V . 

Now consider the event An = An(L) that As occurs for all the strips involved in the 
sausages Si{x,y) for all x,y £ S. Clearly, A N (L) C A N (L') if L < V. 

Definition 3.19 Let 

4 = inf{L : 3L < N and An(L) occurs}, 
with the convention that inf = oo. 

The following result will be proven below. 

Proposition 3.20 

t"> n < C(£ e + l) 2d ( N 2+e + max — !— ^) , (3.35) 

where C > depends only on d. 

This (deterministic) result, together with the following (probabilistic) one yields (3.6), after 
one uses Lemma 3.13 and Borel-Cantelli. 



Lemma 3.21 For all large enough N 

P 





\d+r 




4 




7e 



< 



iV 1 +' 5: 



where c depends only on d and 5 > is independent of N. 
Proof of Lemma 3.21 For L > fixed, we have that 

¥(A C S (L)) < dN max F(x is connected within S by a *-closed path to 5), 

where <S and S are the two larger sides of S. Now the latter probability can be bounded above 
in a standard way by 

1>L 

where A; is the number of distinct *-paths of length I within S and starting at x. This is 
bounded above in a standard way by 7 l and thus 



P(A|(L)) < dNj2( 7N ~ ye ) 1 < cA^ 1_7eL/2 , 



1>L 



for some constant c and all large enough N. 
Then 

F(A C N (L)) < c N 2d+1 -^ L / 2 . 
The result now follows from the observation that {£ e > L} C A^-(L). □ 
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3.5 Proof of Proposition 3.20 

We assume £ e < oo; otherwise, the bound is obvious. We choose the weight function W. For 
b E £>, we make 

w(b) = I 1 : if6ise -g° od ' 

I jj, if b is e-bad. 

We now choose a complete set of paths for S, T. Since £ e < oo, we have that for all x,y E S, 
there will be a (£ e , e)-good path within Se e (x, y) connecting x and y, so we choose one of them 
(according to some arbitrary predetermined order), call it ir x>y , and make 

r = {ir x , y ;x,y E S}. 

We now use the above W and T in (3.34). Let B\ = {b E B : b is e-good} and B 2 = {b E B : b 
is e-bad} =B\Bi. Then 

T «,N < 7^+7^2 + 7^+7^ 2 , (3.36) 

where, for i,j = 1,2, 

3 b ^ ^ W( & ) 

For x,y E S, let Qx, resp. Q y , denote the £ e x 4 square contained in S^(x, y) with x, resp. y, 
as one of its corners. 

Remark 3.22 Notice that for every x,y G S, all bonds ofir X;y \ (Q x U Q y ) are e-good. 

Given b,b' G B, let M(b,V) = #{(x,y) G 5 x S: 6, 6' G 77^}; see Definition 3.17. 
Estimation of r^. 

t£ < A^max^A/XM')- 



beB 

b'eB 



Now for every b,b' E B 



N{b,V) < #{(x,y) E S x S : b,b' E S L (x,y)} 

< #{(#, y) E S x S : a, a E r) x , y for some a,aEB: dist(a, 6) V dist(a', b') < 4} 

< M(a,a'). 

a,a': dist(a,6)Vdist(a',fe')<^E 

where dist is the usual Hausdorff distance between sets. Thus 

t£ < lfN e M N (3.37) 

where M N := A^ d max aeB XL'eB -M( a , a ')- 

To estimate Mjv, we start with the observation that since our paths are described in an 
oriented way, we must specify which of a or a' is traversed first and in which direction. Given 
a = (w, z), we have 

A4(a,a!) = ^^#{(x, y) : a, a' E r\ xy in the order w, z, w', z'} (3.38) 

a'=(w',z')eB a'eB 
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■ a,a' e r] XjV in the order z, w, w' , z'} (3.39) 

a'eB 

+ y) '■ a, a' G r\ xy in the order w', z', w, z) (3.40) 



a'eB 



+ y) '■ a,d & r\ xy in the order vJ , z' z, w}. (3.41) 

a'eB 

We estimate the sum in (3.38). The estimation for the ones in (3.39-3.41) is similar. Let j 
be the coordinate where w, z differ, that is z.i = Wi if % ^ j and Zj = wj ± 1. Then the ordering 
imposes that z\ — w[ — Wi if i < j. The sum in (3.38) can then be decomposed as follows. 

A<'(a,a'), (3.42) 

k=j a'eA k 

where Ai'(a, a') = #{(#, y) : a, a' G r] xy in the order w, z, w', z'} and 

Afc = {(w 1 , z') G B : z\ — w\ — Wi, if i < j; z\ ^ Wi, if j < i < k; z[ — w^, if k < i < d}. 

It is clear that |A fc | < N k ~ j+1 . Now, for a' G A k 

M'(a, a) < #{x e S : Xi = Wi for i > j} x G S : yt = z[ for i < k) 
< N j N d ~ k+1 . 

Thus (3.42) and (3.38) are bounded above by dN 2+d . After a similar reasoning for (3.39-3.41), 
with the same bounds, we finally get from (3.37) that 

t£ < Adl 2d N 2+e . (3.43) 

Estimation of rf 2 . 
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< 7V e - d+1 max ^ A/"(6,6'). 



beB 

b'eB 2 



By Remark 3.22, if b' G $2 is in ir Xty G T, then 6' must be either in Q x or in Q y (see definition 
right above Remark 3.22). Thus 

J\f(b, V) < #{(x,y) <E S x S : a <E r] Xjy for some a G B and dist(a, b) V dist(x, b') < £ e } 
+ y) <E S x S : a G Tfc^ for some a G £> and dist(a, 6) V dist(y, 6') < 4} 

< ^ [J(a,z)+J(a,z)], 

aeB.zeS: dist(a,6)Vdist(z,o')<^ e 

where 

J{a, z) = #{x G S : a G r?^}, j(a, z) = #{y G S : a G 77^}. 
We conclude that 



7^ 2 < £ d N e ~ d+1 



max ^ ,7(0, z)W(z) + max ^ ,7(a, z)W(z) 
aeB zes aeB zes 
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< const e d N e ' d+1 



max J' (a, z) + max > J(a,z) 

aeB ^ r,<zn t-~t 



zeS 



aeB 



zeS 



(3.44) 



since 

W(z) : = #{&' G B : dist(z, 6') < 4} < const £ d . 

We estimate the first max term in (3.44). The other one is treated similarly, with the same 
bound. Let a = (u,v). We decompose J7"(a, z) in J'{a,z) and J~"(a,z), where 

J 1 {a, z) = #{x G S : a G ?7 X)2 , with w traversed before f }, 
J" {a, z) = #{x G S: o £ 77^, with i> traversed before u}. 

We estimate max ag g ^2 zeS J '(a, z). The expression involving J7""(a, z) is treated similarly, with 
the same bound. Let j be the coordinate where u and v differ. Then z must satisfy zi = u iy if 
1 < i < j — 1. We conclude that there are at most iV ' - ^ 1 such z's. For each one, if a G r/ X:Z , 
then x must satisfy Xi — u i: if j + 1 < i < d We conclude that there are at most A" such x's. 
Thus, 

max V J'(a, z) < max JV J ' N d ~ j+1 = N d+1 . 

aeB ^ l<j<d 



zeS 



We conclude that 
Estimation of r£j. 
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t£ < const £ 2d N 2+e . 



< (max^-| TV-**- 1 max VMM' 



(3.45) 



We now estimate the max of the sum above, in much the same way as we estimated 
maxbgg ^ 6 , ei?2 A/"(6, 6') above. By Remark 3.22, if b G £>2 is in tt x ^ G T, then b must be 
either in Q x or in Q y . Thus 

&') < y) <E S x S : a G for some a G £> and dist(a, 6') V dist(:r, 6) < 4} 

+ y) <E S x S : a G 77^ for some a G £> and dist(a, 6') V dist(j/, 6) < 4} 

< ^ [.7(a,z)+^(a,z)]. 



aGB.zGS: dist(a,fe')Vdist(z,6)<4 



Thus, 



max V jV(6, 6') < const £ d 

beB-2 J 
b'eB 



max V J(a, z) W(a) + max V ^(a, z)W(a) 



aeB 



zeS 



aeB 



< const t 



2d 



max > J7"(a, + max > J7"(a, 2) 



zeS 



aeB 



(3.46) 



where 



W(a) := #{&' G B : dist(a, b') < Q < const i d . 
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We estimate the first summand within square brackets in (3.46). The second one can be 
similarly estimated with the same resulting bound. 

max ^ J(a, z) < max ^ J' {a, z) + max ^ J" {a, z) (3.47) 

2 aeB Z aeB % aeB 

and we estimate the first summand within square brackets in (3.47) only. The second one can 
be similarly treated with the same bound. Let z G S be fixed and j be the coordinate where 
u and v differ, where (u,v) = a. Then u must satisfy Ui = Zi, if 1 < % < j — 1. We conclude 
that there are at most N d ~i +1 such w's. For each one, if a G r] XjZ , then x must satisfy Xi = Ui, 
if J ' + 1 < i < d. We conclude that there are at most N j such x's. We then conclude that 

d 

max J'( a , z ) < max Yl J '( a > z ^- d N<l ~ 3+l N * = d Nd+1 > 

Z a£B Z j=l a=(u,v)<EB 

u and v differ in j 

which eventually yields 

7& < const i\ A (max-^r) • (3.48) 



x£S Uj(x) J 

Estimation of t% 2 . 

t" < (max4rl N ~ d max V Nib, V). (3.49) 

22 ~ V xdS Ujix) I bdB 2 ^ 

By Remark 3.22, if b, b' G £> 2 is in n x>y G T, then 6' must be either in Q x or in Q y (see definition 
right above Remark 3.22). Thus for b G B 2) we have 

A/"(M') < E e Qx}+Hb,bf G Q y } + l{b g G Q„}+1{6 G Q v , f G Q x }]. 

6'ei3 2 x,y£Sb'£B 2 

Now 

£ ]T 1{6, b'eQ x }<J2J2 ^ b e Q*t E ^ G Q*}- ( 3 - 5 °) 

x,y&Sb'£B 2 yes xeS b'eB 2 

The two inner summands in the left hand side of (3.50) are uniformly bounded by const if , so the 
left hand side of (3.50) is bounded by const Pf N d . For similar reasons, the same bound holds 
for E^esE&'^UM' e Q„}, E^s^e&l^ ^ Q^fc' G Q„} and E^esE^s, e 
Qy,^ G Qa;}, and thus, from (3.49) 

r 2 % < const ^ d ( max-)— J . (3.51) 

22 ~ e \ xes lo(x) J v ' 

The result of Proposition 3.20 now follows from (3.43), (3.45), (3.48), (3.51) and (3.36). □ 



4 Decay of the annealed return probability for random 
walks on Z rf 

We go back to the study of Markov chains taking their values in Z d . Let uo : Z d — > and 
define the Markov generator 
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where the sum is over sites y which are nearest neighbors to x. 

As in Section 2, {X t , t G R+} will be the coordinate process on path space (Z d ) R + and 
we use the notation P£ to denote the unique probability measure on path space under which 
{X t , t G M + } is the Markov process generated by (4.1) and satisfying X = x. 

As in Section 3, we choose the family {u(x), x G Z rf } at random, according to a law Q on 
(R;f d such that 

the random variables {ui(x), x G Z d } are i.i.d. ; 
uj{x) < 1 for all x; 
Q(w(0) < a) ~ a 7 as a | 0, (4.2) 

where 7 > is a parameter. 

Remark 4.1 We note that this generator has the same form asQ w in (2.1) by making u(x,y) = 
uj(x) /\uo(y), and also the same form as C U,N in (3.1), but in infinite volume. There would also 
be similar results for uj defined on edges, instead of points, with i.i.d. values for different edges, 
and the same technique would apply. 

Remark 4.2 lfu(0) were a Bernoulli random variable, then we would have a random walk on 
a (independent, site) percolation cluster (provided we started in an infinite occupied cluster). 
See [6]. 

In the sequel Q.P£ will be used as a short hand notation for the annealed law defined by 
Q.P£[ •] — J Px \ ' ]d<Q.(ou). We are interested in estimating the decay of the return probability 
under Q.P W , Q.P w [X t = 0], as t tends to +00. It is actually quite easy to derive lower bounds 
for Q.F^lXt = 0]. Indeed, on one hand, one can use the comparison lemma 2.2 with the usual 
nearest neighbor random walk on Z d to prove that 

Q.V^lXt = 0] > ct~ d/2 , (4.3) 

for some contant c that depends on the dimension d. There is another way to prove (4.3), 
as follows. It is known [3] that, under Q.P W , X t satisfies the central limit theorem. Together 
with the reversibility and the translation invariance of the law Q, the C.L.T. implies (4.3) (See 
Appendix D, in [6]). 

On the other hand, for any realization of u, the first jump of X t follows an exponential law 
of parameter ^ o w(0) Aw(?/) < 2duj(0). Therefore 

Ki x t = 0] > ^[X s = 0,Vs < t] = e -^v~o^W > e -2<M0)t_ 

Taking expectation w.r.t. Q and using the condition (4.2) on the law of a>(0), a simple compu- 
tation leads to a lower bound of the form 

Q.P W [X t = 0] > cr 7 . (4.4) 

As is indicated in the next statement, these lower bounds turn out to be of the correct 
logarithmic order. 
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Theorem 4.3 



lim 



logQ.P"[X, = 0] 
log* 




(4.5) 



Remark 4.4 From the point of view of statistical mechanics — here the statistical mechanics 
of a disordered system — we consider Theorem 4-3 as a (nice) example of a dynamical phase 
transition. 

Remark 4.5 Such tools as Sobolev embeddings, isoperimetric or Nash inequalities of constant 
use for estimating transition probabilities of Markov chains, see [2], cannot be directly applied 
here because of the lack of ellipticity of the transition rates uj. Thus (4-5) is also an example of 
exotic 'heat kernel decay' for a non uniformly elliptic generator. 

Remark 4.6 A fruitful technique to handle r.w.r.e. is to isolate the effect of the fluctuations of 
the environment uj in a given scale. See for instance random walks in Poisson environments [11] 
where one single eigenvalue dominates the rest of the spectrum. There does not seem to exist 
such a separating scale in our model. 

In view of (4.3) and (4.4), only the upper bound is missing in the proof of (4.5). 
We use spectral theory. We rely on a trace formula similar to the one obtained in Section 2 
and on our spectral gap estimates from Proposition 3.20. 

4.1 Trace formula 

We express the annealed return probability as a trace. The argument is the same as in Section 2, 
except that we restrict ourselves to computing the trace on cubes whose radius can be chosen 
as a function of time. This is possible because rates are assumed to be uniformly bounded. 

Let £ > 0. In the sequel, we shall use the notation N = t^ 1+ ^/ 2 . (In fact, N should be 
defined as the integer part of *( 1+ ^/ 2 , but, for notational ease, we will omit integer parts.) 

Let B N = [-N, N] d , be the box centered at the origin and of radius N. Let C U ' N be the 
restriction of the operator iy to B N . Thus C U ' N is defined by 



where the sum is now restricted to neighboring points x and y in and we impose peri- 
odic boundary conditions. —C W,N is then a symmetric operator. We denote by {\f(BN),i e 
[1, #B N }} the set of its eigenvalues in increasing order. 
Let tn be the exit time of X t outside B N . 

We compute Q.F%[X t = 0] using the translation invariance of the probability Q. Since 
Q.P£[X t = x] does not depend on x, we have 




(4.6) 




1 



X] t < T 2N ] + 



1 



£ Q.^ x [X t = x;t>r 2N ] 



xeB N 



#B N 



xeB N 
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^ ^ xeB 2 JV ^ N x&B N 

If under P£, x G -Bat, we have t > t 2 n, then the process must have left the ball x + before 
time t. Since the probability Q.P£[3s < t s.t. X s x + Sjv] does not depend on x, we have 
that Q.P£[t > r 2JV ] < Q.Pft[t > 7*]. 

We note that XLeB 2 iv ^xt-^* = x;t < r 2N ] is the trace of the semi-group of the process X t 
killed when leaving the box B 2 n, i-e., with Dirichlet boundary conditions outside B 2 n- It is 
therefore dominated by the trace of exp(t£ w ' 2Ar ), that is 

Y, K[Xt = x-t<r 2N ]<Y,e- X " {B ™ )t . 

x<=B 2 n i 

Thus, we have proved that 

Q.F%[X t = 0] < -i- VQ[ e - A ^ B ^*] + Q.F"[t > 7*]. 

From the Carne-Varopoulos inequality, it follows that 

P"[t > rjy] < 2tN d - 1 e-^ + e _c ', (4.7) 

where c is a numerical constant, see Appendix C in [6]. With our choice of iV = t( 1+ ^/ 2 , we get 
that P w [t > r^] decays faster than any polynomial as t tends to +oo. 
Thus Theorem 4.3 will be proved if we can check that 

limlimsup 1O g^ e " W U0V^- 7 V (4.8) 
5 ->o t _> + J log* \2 7 

4.2 Min-Max 

C is a constant that depends only on d and Q. For constants depending on other parameters, 
we indicate it. 

Let us first recall the lower bound on the first non trivial eigenvalue of an operator of the 
form C W,N . In Section 3, we proved that 

- * r <c( N 2+£ + sup -L-} (4.9) 

In (4.9), e is any positive number; C is a constant depending on the dimension only; d^ is a 
measure of the set {x e B N : uj(x) < N~ £ }. 

With the notation of Section 3, Proposition 3.20, d® = (4 + l) 2d - (But note that l e depends 
on JV.) Thus d^ is a random variable, i.e., depends on e, N and also u>. 

Using the properties of Q, we get that, for some constant c, that depends on Q only, we 
have 

> A) <cN-^, (4.10) 
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where A can be chosen such that A > — and N is supposed to be large enough. (How large 
depends on the dimension only.) A proof of (4.10) can be found in the proof of Lemma 3.21. 
From the min-max caracterization of the eigenvalues of symmetric operators, we have 



^iVi(-Sjv) = max min - 



lE^ fljf ki)A W (y)] [f(x)-f(y)} 



2 



A--./* / 2 V,.. /;x P(.r) 

where the 'max' is computed on choices of % functions defined on and the 'min' is computed 
on functions / such that, for all j G Y1x&b n f( x )fj( x ) = 0- 

Thus, in the computation of \^ +1 (Bn), we may impose at most i different linear constraints 
on the test function /. We consider two kind of conditions. 

Let k G N*. We chop Z d into a disjoint union of boxes of radius k, say Z d = U ze %d~B z , where 
B 2 = (2k + l)z + Bk- We now choose for some of the function f/s, the indicator function of 
the boxes B 2 that intersect Bn, i.e., we require that 

E /(*) = °> 

x£B N nB z 



for all z G Z d such that B N fl B 2 ^ 0. The number of such z's is at most 

/ _! . \ -+- ! - IK - ! ^ " 

n 2 = 



Clearly, 



Ik 



E/ 2 (-) = E E m 

x£Bpf z x£B N nBz 



and 



E c(x)Ao;(y)(/(x)-/(y)) 2 >E E K*) A c,(y)] [/(x) - /(y)] 2 . 

x^y&B^f z x^y^BpfPiBz 

Therefore 

E^ % kWA^)][/W-/(y)] 2 . £^ 6fl W w (*) Mv)l [/(*)" /(v)] S 

> mm ■ 



SxeB]v f 2 ( x ) z ^x<=B N nB z 

where, for each z G Z d , J2xeB N n& z 

f(x) = 0. 

Next, let us choose n x points in B N , say Si,...,S ni . We choose for some of the f/s, the 
indicator function of the points Sj and their neighbors in B N , i.e., we specify that f(x) = 0, 
for x G {Si, ...,S ni } or x ~ <5j, for some j. This recipe leads to, at most, (2<i + l)rii different 
conditions. We note that, for such a function /, the value of the Dirichlet form 



E [u,(x)Au,(y)}[f(x)-f(y)} 

x~y&B N 



2 



does not depend on the value of ou(Sj) anymore. Therefore, we may assume that oo(Sj) = 1, for 
3 e [l,ni]. 
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Thus we see that, if i > n 2 + (2d + l)ni, then 

K+i(B N ) >minAf(S^nB z ), 

z 

where u is a new environment obtained by modifying the value of a; to 1 on all points Si and z 
ranges through those points in Z d such that B 2 intersects B N . 

We now choose for Sj the points in B N where u achieves its lowest values. Let us use (4.9) 
to estimate each eigenvalue \%(B N fl B 2 ): 

TZT^T < C (k^ + sup -L) d-. (4.11) 

We used ci^ as a uniform upper bound for the minimal side length of strips for which the 
event A^(L) in Definition 3.19 occurs. 

sup™g Bjv l/cu(x) denotes the maximal value of l/u(x), i.e., 

sup . . = max{/i : #{x G Sjv : u>(x) — 1/h} > n\ + 1}. 



4.3 Proof of Theorem 4.3 



Remember that we have already chosen some parameter £ > (that we want to choose close to 
and which is related to t by N = t^ 1+ ^/ 2 ), and another parameter e > which is arbitrarily 
close to 0. We need a third parameter a G (0, 1). The constant A in (4.10) is at our disposal. 
We also still have to choose n\ and n 2 , depending on % and such that % > n 2 + (2d + l)ni. 
Write 



< 



+ £ 

i 

< (2N + l) d Q[d? >A] + J2 



e -V(B N )t. x«(B N ) > N~ 



e -X?(B N )t. d N < ^ A . (jBiv) < N - e ( ^ 



df < A and A^-Bjv) < N~ 



•a/d 

w 



Using (4.10), we see that we can choose A in such a way that 

log[(2iV + l) d Q(d^ > A)] ^ nxi (d 

hm sup : < V 7 

t-^+oo \ogt \2 



•a/d 

W 



(4.12) 



(4.13) 
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An easy computation shows that 



lim sup 



-N~ 



■a/d 



logt 



< 



d_ 

2a 



S 2 2 



(4.14) 



Let us now bound the last term in (4.12). Assume that < A and X^(B N ) < N 



From (4.11), we must have 



■ f a/d 



N e+2 i~^ < C k 2+£ + 



ni 1 

sup —— 

x£B N UJ(X) 



A. 



We choose i%2 = i a and assume that i is large enough, how large depending on the dimension, 
a and 7 only, which we may do. Then k 2+£ < N 2+£ i s . Therefore, we must have 

A^ +2 2 d < C sup — -, 

with a possibly different value for C. 

From now on, we deal separately with the cases of large or small values of 7. 

Case 7 > tj. We then choose e < 2- and a — 1 - ^ + §. 

The computation goes as follows (the value of C changes from line to line) 



s E 1 

i 

< CQ 



C < A and A^(Bjv) < iV" 



ja/d 



< 



E' 



ni 1 

C sup — — — 



> iV £ + 2 



_ 2a 



ni 

C sup 

xeB N 



uj{x) 



> N' 



£+2-2a 



# \ i : C sup 



o;(x) 



> N' 



e+2-2a 



C sup — - 



> N- 



e+2-2a 



= C(2N+1) C 



uj(x) 



> N' 



e+2-2a 



< CN d N~^ £+2 ' 2a) = CN-'^r, 



(4.15) 



where the second inequality follows because % < (2N + l) d , and the third one because n\ 
(i - i a )/{2d+ 1) and a < 1. Thus we deduce from (4.15), (4.12), (4.13) and (4.14) that 



logQfE^^] ^ d 

lim sup < — 

logt 2a 



e e £ 

£ H h — 

s 2 2 



Let e tend to and then £ tend to to deduce (4.8). This ends the proof of Theorem 4.3 in 

the case 7 > d 



Case 7 < |. Let 5 G (0, 7), to be chosen later. We have 



rff < A and A?(fl*) < AT- (^) 

< tf(2a-e-2)8 £^ ,j- 



<E.Q 
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since % < (2N + l) d . Remember that n 2 = i a is much smaller than n\ — (i — i a )/(2d + 1) for 
large values of i, say i/(4d + 2) < n\ < i/(2d + 1). Let x , ...Xj, ...X( 2 N+i) d -i be an enumeration 
of the points in B N such that the sequence co(xj) is increasing. Thus 



ml 1^1 
sup — — = — < — v 

xeB N UJ{X) UJ{X ni ) CO {Xi/(4 d+2 ) ) 



Therefore 

i 

where c§ = 



m 1 
sup — - 

X& B N 0J[X) 



<{Ad + 2) 1 

xGBm 



Uj(x) 



c 5 (4d + 2)(2N + l) d , 



i 

ui(x) 



. Note that Q 



is finite and does not depend on x. Therefore 



< A and \?(B N ) < N~< ( ^ 



< C(5 A^ (2a ^- 2)5 (2iV+ l) d . 



Gathering this last inequality with (4.13) and (4.14), we get that 

logQ^e-W")*] 



lim sup 



log* 



d 



d 



< max ^0-7;^)^ + ^ + 



e , e£\ 1+e 



[d + (2a - e - 2)5] \ , 



(4.16) 



with a G (0, 1) and 5 G (0, 7). First replace 5 by 7. Then let e tend to and choose a = ^z§^- 

The upper bound in 
and conclude that 



The upper bound in (4.16) becomes max fj — 7; : ^(d — 27) + ( - 1 ^^ 7 . Finally let £ tend to 



,. ,. logQ^e-W")*] d 

lim lim sup < 7, 

€-0 t^+oo log* ~ 2 



and Theorem 4.3 is now proved in the case 7 < f . □ 



5 Quenched decay of the return probability 

In this section, we investigate the quenched decay of the return probability. Model and notation 
are the same as in Section 4: a random walk among i.i.d. random conductancies with a power law 
with an exponent 7. Now we are rather interested in the asymptotics of the return probability 
Pq [X t = 0] in Q probability. Let us set a c to be the best exponent a such that 

Q\P%[X t = 0] < r Q ] -> 1 as t -> 00 . (5.1) 

From Theorem 4.3, it is clear that a c > | A 7. We can do better in the case 7 < | : 
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Theorem 5.1 For any 7 < f then a c > 7. 



Remark 5.2 Although rather unsatisfactory — because it does not give the true vaalue of a c 
- Theorem 5.1 shows that the typical decay of the return probability is strictly faster than the 
averaged decay. Such a situation is sometimes called in the litterature a 'high disorder regime '. 

Remark 5.3 The proof of Theorem 5.1 actually yields the lower bound 



There is no reason to believe that this bound is sharp for a given value of 7. Notice however 

|, we get the inequality a c > \, 



that, in the regime 7 — > | , we get the inequality a c >\, which seems to be sharp. 



Let us sketch the proof: we use the fact that, with large Q probability, the origin lies in an 
infinite percolation cluster, say C, of 'good' sites, where u is bounded from below. Estimates 
on the return probability for random walks on percolation clusters have been proved in [6] (See 
also [1] ). One strategy would then be to try to couple the random walk in the environment uj 
with the random walk on C: we have no idea on how to do that. We rather rely on spectral 
theory to compare the behaviours of the eigenvectors for the two random walks. Note that from 
the results of [6] follow precise estimates on the eigenvalues of the discrete Laplace operator on 
C. The core of the proof is to show that eigenvectors of the generator of the random walk in the 
environment u, when they correspond to small enough eigenvalues, are concentrated outside C, 
and therefore do not contribute too much to the asymptotics of the return probability as soon 
as the random walk starts outside C. 



5.1 Proof of Theorem 5.1. Step 1 

Let a < I ■ Choose two parameters e > and £ > 0. We shall use the notation N = 
i^+O/ 2 . (j n f ac ^ ]\[ should be defined as the integer part of t( 1+ ^/ 2 , but, for notational ease, 
we will omit integer parts.) All the limits to be taken are to be understood as t — > 00 or, 
equivalent ly iV — > 00. 

Let C w be the largest connected component of the set {x E Z d : uj{x) > N~ £ }. We assume 
that iV is large enough so that Q[u(x) > N~ £ ] becomes larger than the critical percolation 
probability on Z d . Then C w is the unique infinite connected component of the set {x E Z d : 
uo(x) > N~ £ }, see [5]. We denote by C% the largest connected component of the intersection 
C u n B N , where B N = [-N, N] d . 

In the next step of the proof, we will define a set of environments, denoted Q^, such that 
®p N ] -> 1- We further have the property Q[f§|] ->• 1. 

Calling {\^(B N ),i e [1,#Sjv]} the eigenvalues of —C U > N in increasing order, and {ip?,i G 
[1,#Bn]} the corresponding eigenvectors with due normalization in L 2 (B^), a very similar 
computation as in Subsection 4.1 leads to the following series of inequalities. 

We first use the invariance by translation of Q. 

QK\x t = 0] > r a ] = Q[ir x [x t = x}> r a \ 
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holds for any x G B^- Therefore 



%[x t = o] > r] = ^-E Q [ p ^ x * = ^ ra l- 

Note that 

*Z[X t =x}> r a ] < QK[X ( = x; * < r 2N ] > t-72] + Q[P£[t > r 2N ] > t~ a /2] , 



where t 2 n is the exit time of B 2 n- Since sup w sup 3, P£ [t > t 2 n] decays faster than any polyno- 
mial, see (4.7), we have 

limsu P Q[P£[X t = 0] >r a ] <limsup— j— V Q[^[X t = x; t < r 2N \ > t~ a /2\. 



xGB N 



We now restrict our attention to those environments belonging to Q N and to the points 
x G Cfj\ 



< 



[X t = x;t<r 2N )>r a /2] 

[X t = x;t< r 2N ] > r a /2; x G C%; £l N ] + Q[Q C N ] + Q[x <£ C%}. 



Since Q[pjv] — > 0, we therefore get that 

limsupQ[P£[X t = 0] > t~ a ] 
< limsup — 4- V QKi x t = x;t< t 2N ] > t~ a /2; x G C%; Q N ] + limsupi 



#5 



iV 



But since 



\*£1 



— > 1 (see step 2 below), we have 
limsupQ[P£[X t = 0] > r a ] 



< limsup — [— V Q[P£[X t = x- 1 < r 2N ] > t~ a /2; xeC%;Q 



N ■ 



From the Markov inequality, we deduce that 

QK[X ( = x- 1 < r 2N ] > r a /2; x G C%; Q N ] < 2t a Q[F" x [X t = x; t < r 2N \; x G C^; Q N ], 
and thus 

limsupQ[P£[X t = 0] > r a ] < 21imsup-^-Q[ V P£[X t = x;t < t 2N ];Q n }. 

ii T — * 



A 1 



Finally we express the probability P£ [X t — x; t < t 2N ] in the spectral decomposition as 



[X t = x;t< r 2N ] = ^-J2 e ' Xf(B2N)t ^^ 2 ' 
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and get that 



t a 

limsupQ[P£[X t = 0] > r a ] < 21imsup — — d 



-Xf(B 2N )t_ 



N 



£(C(*)) 2 ;Ov 



xdCi 



.(5.3) 



Let us pause a little to look at (5.3). It is true that -^g- ^2 xeC ^(ipf(x)) 2 < 1; but if we 

would use this upper bound, we would be left with Q[^2 i e~ x i^ B2N ^ t ], and the best value for 
a would then be 7, as the results of Section 4 show. We have to find a better way. Note 
that terms corresponding to large values of i, and thus large values of \^{B2n), can be easily 
controlled. Thus the main point is to show that ^jj^ X^eC^^^)) 2 ls small enough for small 



i, i.e. we have to prove that eigenvectors corresponding to small eigenvalues are concentrated 
outside Qn. And in fact one would expect this to be true since small eigenvalues arise because 
of small values of u, and these precisely sit outside C%. 



5.2 Step 2. Definition of £In 

The set Qn is defined by two requirements: we ask that for any w G Ojv we have 

(i) e C%. 

The second requirement deals with the behaviour of the random walk on C%: let G 
[1, be the eigenvalues of the discrete Laplace operator on C% as defined in [6]. We will 

also use the notation {4>i,i G [1,#C^]) for the corresponding eigenvectors. We assume that 
the eigenvalues are in increasing order and the eigenvectors are normalized in L 2 (C^) for the 
counting measure. Of course the /XjS and 0jS depend on u and N . 

Let 



7] = and j = N d ~ v . 



2 + d 



Note that since 7 < d/2, then rj < d. We then require that, on £l N , 

■2/d 

(ii) Hj > 



N 2 (\ogN) 8 ( d -v)/ d ' 

The definition of fl^ is now complete and all that remains to be done is to check that 

Q(^7V) -> I- 

That Q((i) holds) = Q(0 G C%) -»• 1 is obvious. 

As for condition (ii), we rely on the results of [6]. Calling P^[Xf = y] the transition 

probabilities for the random walk on C%, we quote from formula (6) of [6]: Q-a.s. on the set 
where C w is infinite 



sup 



1 



A? 



2d 



< c 



(logAQ 

d , jloglogJV 5 
c2 + ° log JV 



where C is a dimension dependent constant, s is arbitrary, and iV > N (u) is large enough. 
(In [6], formula (6) is deduced from the isoperimetric inequality (4), (4) is a consequence of (21), 
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and (21) is proved for both site and bond percolation models with parameter p close enough to 
1, which is our case here. Besides, we replaced e(N) by its value e(N) = d + 2d lo f^f^ N , noticing 

that (4e(iV)//3 2 ) e(iV)/2 then behaves like a constant.) 

We then choose x — y, sum over x G C%, and express the result as a trace to get that 



2d 



(logAQ 

d i j log log N 
S2 t0 logJV 



e-^ s < 1 + #C^C 

i 

Therefore 

,ye b i + ^fL-iv d . , iogiogiv 

g 2 +a log JV 

Take now s = -gjQog N)*^/* . Then 



2d 



\ r jtc» (logA ° o 

■ ^W^N d , ^ loglogAf ^ U ) 
J S 2+ a logJV 



so that e ^ s — > 0, and we have proved that Q-a.s. on the set where is infinite, for large 
enough N, condition (ii) is fullfilled. 

Finally we already used the fact that Q[^g^-] — > 1 that should be justified: from the result 
of Appendix B of [6], we know that the expected density in B N of the component of C u fl B N 
that contains the origin goes to 1 as N — > oo, and Lemma 3.12 implies that the expected 
density in of the largest component of C w D B^ goes to 1 as N — > oo. Thus the component 
of fl Bn that contains the origin and C% coincide for large N and its density tends to 1. 



5.3 Step 3. Spectral analysis 

Assume that uo G fijv- 

We bound the term XLe^^f ( x )) 2 i n (5-3) in two steps by writing that 

where P J is the projection on the subspace of L 2 (C%) spanned by the eigenvectors ((pi, % G [1, j]). 
On one hand, since XLe£iv( < M ;r )) 2 — 1> then 

i ^ N xec% 

s E#^EWM) 2 

- E#^E(«*)> 2 
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On the other hand, for any function / on C%, we have 



E (f(*)-f(y)) 2 , 



this last expression being the Dirichlet form of the random walk on C%. Since u(x) > N £ on 
C%, we get 

E U'( x )-f(y)) 2 < N£ E M*)Ac,(y))(/(x)-/(y)) 2 , 

this last expression being now the Dirichlet form of the random walk on B 2 n- Since ipf is an 
eigenvector, 

— 1— E (^(^)Aa;(y))(C(^)-C(y)) 2 = Ar(S 2J v). 



So 



#5 



From these two estimates, we deduce that 

. _ 1 . _ ATE 



(5.4) 



5.4 Step 4 

Combining (5.3) and (5.4), we see that Theorem (5.1) will be proved once we have checked that 
-g-N^ -> and that 



t a 



-N £ Q 



N 



0. 



(5.5) 



We recall that a < f N = t^l\ V = and N > (^) 2 (log A0~ 8(d_r?)/d on Q N . 

It is then immediate to see that -J^-N d ~ r ' — > 0. Besides, (5.5) will hold for any a < - 1+ ' 



and some e > if 



hmlimsup < 1 — 7. 

v log* ~ 2 



2 l+d/2 



(5.6) 



But using the inequality A^e Al * < ye 2 A ** 5 we get 



lim lim sup 



Q[E,Ar(j? 2j y)e- A "^) f ] 
log* 
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< 



1 + lim lim sup 



Q[J2 l ^ Xf(B2N)t/2 ] 



log* 



< 




by (4.5). □ 
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